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Abstract - We prove the equivalence between the hard-sphere Bose gas and a system with 
momentum-dependent zero-range interactions in one spatial dimension, which we call extended 
hard-sphere Bose gas. The two-body interaction in the latter model has the advantage of being a 
regular pseudopotential. The most immediate consequence is the existence of its Fourier transform, 
permitting the formulation of the problem in momentum space, not possible with the original 
hard-core interaction. In addition, in the extended system, interactions are defined in terms of 
the scattering length, positive or negative, identified with the hard-sphere diameter only when it 
is positive. We are then able to obtain, directly in the thermodynamic limit, the ground state 
energy of the strongly repulsive Lieb-Liniger gas and, more importantly, the energy of the lowest- 
lying super Tonks-Girardeau gas state with finite, strongly attractive interactions, in perturbation 
theory from the novel extended hard-sphere Bose gas. 
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Introduction. — Low-dimensional systems have been 
very popular among theorists for a number of decades 
now [HH], not without a reason. On one hand, many 
of these allow for an exact treatment. The Bethc ansatz 
solves a number of one-dimensional problems exactly, 
such as the Lieb-Liniger gas [3] or Heisenberg's model 
[3], while systems with supersymmetric Hamiltonians [S], 
such as Sutherland's model or the attractive Lieb- 
Liniger model [7], allow for a trivial evaluation of their 
ground state wave functions [T]. On the other hand, one- 
dimensional systems have physical properties of interest, 
and can become strongly correlated, as is the case of the 
Tonks-Girardeau gas [51I2 - a one-dimensional system of 
impenetrable bosons sharing many of its properties with 
the ideal Fermi gas - which was experimentally realized 
in [TO]. More recently, the metastability of the so-called 
super Tonks-Girardeau gas - a strongly attractive one- 
dimensional Bose gas with no bound clusters ~ has been 
proposed |llfll3| and subsequently experimentally realized 
[14] . This unique system attracts much of current interest 

[TOHIl]- 

There is strong theoretical evidence [H] on the close 
relation between the lowest-lying super Tonks-Girardeau 
(sTG) gas state and a system of onc-dimensional (ID) 



hard-spheres in the low-density limit. This is indeed very 
appealing, since the hard-sphere gas is extremely simple 
to solve with a variety of methods (see e.g. |8]), while a 
good description of the ID Bose gas in the sTG regime re- 
quires, in principle, accurate numerical calculations, such 
as the diffusion and variational Monte Carlo employed in 
[13]. This evidence was then used to conjecture that, in 
a harmonic trap, the sTG and ID hard-sphere Bose gases 
are equivalent [19]; we call this statement Astrakharchik- 
Girardeau conjecture for the trapped case. 

In this Letter, we prove a theorem which states the 
equivalence between ID Bose gases interacting via hard- 
sphere potentials and certain momentum-dependent con- 
tact interactions - which we also define in general - in 
the untrapped case. Our theorem represents a weaker, 
though exact version of the Astrakharchik- Girardeau con- 
jecture in free space. The two-body interactions in our 
new model, which we call extended hard-sphere (eHS) sys- 
tem, are Fourier-transformable and we are therefore able 
to write down the system's Hamiltonian in momentum 
representation. For the Lieb-Liniger gas, we obtain the 
energy of both the ground state with strong repulsive in- 
teractions and the lowest sTG state with strong attraction 
to first-order in perturbation theory, departing from the 
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exactly-solvable eHS model as the zero-th order reference 
Hamiltonian. Last but not least, we obtain the so-called 
Tan relations j20H28] for the Lieb-Liniger gas and use them 
to calculate Tan's contact - related to short-distance cor- 
relations - in the homogeneous and trapped cases, the 
latter within the local density approximation. 

System Hamiltonians. — The Hamiltonians for the 
hard-sphere (HS) and Lieb-Liniger (LL) systems with N 
identical bosons of mass m have the form 



H 



N 



/ — + y 

i—l i<j—l 



V{xi - Xj). 



(1) 



For the LL model, V{x) = ghL^ix), 5ll = —2h^/ma con- 
stant, and a the two-body scattering length; for the HS 
model, V{x) = (oo) for |a:| > a {< a), with a > a 
constant hard-sphere diameter. 

Two-body problem. — We begin by considering the 
two-boson problem. After separation of ccnter-of-mass 
{X = {xi +X2)/2) and relative (x ^ xi — X2) coordinates, 
Hamiltonian ([T]) reads 



H 



(2) 



with ^ = m/2 the reduced mass of the two-boson system. 
The stationary Schrodinger equation Hij: = Etp at positive 
energies E = h'^k'^/2fj. is solved by t/j{x) = sin(fc|a;| -I- 9a)- 
The wave functions of the HS model are given by ip^x) for 
|a;| > a and are zero for \x\ < a, while for the LL model 
they are given by ipix) for all x. The phase shifts, with 
self-explanatory subscripts, are given by 



tant^LL = 

tan^HS ~ — tan(fca). 



(3) 
(4) 



Note that if, in Eq. ([3]), jll is made momentum- 
dependent, ^LL ^ 9{k), then the HS phase-shifts (|4]) are 
recovered by choosing 



g(k) ~ cot[ka). 



(5) 



The above relation, carefully stated, provides the desired 
mapping between HS and Dirac delta interactions in one 
dimension. 

Definition of momentum-dependent interac- 
tions. — Any analytic function of an operator is defined 
by its expansion in powers of the operator |29j . Given that, 
define an even analytic function of the momentum oper- 
ator f{k), with k = p/h. The action of f{k) on a plane- 
wave is given by /(fc)e*''^ = f{q)e'^''^. It is also necessary 
to consider its action on functions of the form sin((7|x|), 
with discontinuous derivatives at the origin, which is 
given by Er=o(-1)"/('"HO)(9.)'" sin(g|x|)/(27i)!. Ah 



even derivatives of sin(g|a::|) include undesired Dirac deltas 
S{x). A simple way to deal with this problem is to restrict 
the action of / to positions a; > or x < 0, which avoids 
Dirac deltas and has no influence on differentiablc func- 
tions cos{kx)) at the origin. We can therefore define 
momentum-dependent contact interactions as follows: 

Definition. Let / be an even, analytic function. A 
momentum-dependent contact interaction is an operator 
Wk with the action 



[Wk^]{x) ^ 5{x) 



lim / 



dx 



(6) 



With the above definition, Wk is not Hermitian. This is 
a general property of momentum-dependent pseudopoten- 
tials, shared by the famous partial-wave pseudopotentials 
of Huang and Yang [3^ . 

Equivalence between hard-sphere and contact 
interactions. — From the above considerations, the 
proof of the equivalence between the HS states and the 
states corresponding to momentum-dependent contact in- 
teractions for the two-body case is immediate: 

Theorem. Let H be defined by Eq. with 



V{x) = 5{x)g{k), 



(7) 



where g{k) = — fPk cot ak, and the limit x — > 0^ as in Eq. 
^ is assumed. Let ipk : M C be the bosonic scattering 
wave functions of H, at energies fi,^fc^/2/i, k real. Then, 
the restrictions (j)k oi ipk to V = M. — {—a, a), a > 0, are 
the bosonic scattering wave functions at the same energies 
for the problem of two hard spheres of diameter a. 

Remarks: (i) a > is necessary for the mapping, al- 
though Hamiltonian ^ with potential ([7]) is well-defined 
for a < 0, so we regard a as the scattering length of 
the model, which we call extended hard-sphere (eHS) 
model; (ii) it is fundamental in the above theorem that 
the energies are positive, since there exists an unphysi- 
cal, infinitely-bound state which becomes the identically 
zero function in the restricted domain of the HS wave- 
functions. 

Many-body problem. — Fortunately, the many- 
boson problem with contact interactions is exactly solv- 
able by means of the Bethe ansatz. This leaves us with 
the equivalence of the HS and the extended model, stated 
as follows 

Theorem. Let H be defined by Eq. ([T]), with 



where fc, 



V{xi - Xj) 



S{xi 



(8) 



-h^k cot ak, and the 



fc,)/2,^g(fc) 

limit Xi — Xj — > 0+ as in Eq. ([6]) is assumed. Let 
'0fci,fc2,....fcN ■ — !• C be a bosonic scattering wave 
function of H with energy E = X^iLi fi^kf /2m and k^ 
real, i = l,...,iV. Then, the restrictions (l>ki,k2----,kN of 
'0fci,fc2,....fejv to the domain T> = 1 

= {{xi,X2,...,xn) 6 



1>N\ 



with 



i<3 = 



<«}, (9) 



p-2 



Exact equivalence between one-dimensional hard-sphere and zero-range potentials 



are the bosonic scattering states for the A^-boson problem 
with HS interactions of diameter a > at the same ener- 
gies. The same conclusion holds valid in a box of length 
L > Na. 

Proof. Showing that Bethe ansatz (BA) wave functions 
constitute the bosonic scattering states of Hamiltonian ([1]) 
with interactions ([5]) is trivial: the momentum-dependent 
interactions only multiply each term of the BA by a con- 
stant (which depends on the relative momenta, obviously) . 
Since the potentials have a zero range, scattering occurs 
without diffraction [31] and therefore the model is exactly 
solvable via BA. For the finite box case, the same holds ev- 
idently true. Now, because the two-body scattering phase 
shifts for the eHS model are identical to the HS phase 
shifts, the BA equations are identical, too. In the HS 
case, the saturation density, for which the ground state 
energy at finite densities (in a finite box) diverges, is given 
by N/L = 1/a. At that point the wave functions simply 
do not exist. Now assume that N/L < a and take the 
restriction (j)ki,...,kN of a BA wave function V'fci,. to 
v. Since the BA equations for the two models are identi- 
cal, 4>ki,...,kN s-re the scattering wave functions for the HS 
problem. QED. 

Momentum space. — We now turn to Fourier- 
transform any general momentum-dependent contact in- 
teractions, including the particular eHS interactions. Our 
starting point is the interaction in terms of the bosonic 
field operators -0, 

Wk = ]^j i^Hx)i}\x')Wk{x-x')i^{x)i^{x')dx'dx, (10) 

The field operators are expanded in the plane-wave basis 
as ^{x) = J2p flps'^^, with L the length of the system, 
and flp the bosonic annihilation operators in momentum 
space. Inserting this expansion into Eq. (|10p . we obtain 
the momentum representation of Wk^ 



Wk = 



- E 

pi,P2,q 



(11) 



where pi^2 = (pi— P2)/2. In particular, Eq. (fTTj) applied to 
the extended HS interactions provides, for the first time, 
a simple momentum representation for a very singular po- 
tential. 

Strongly-interacting delta Bose gas as a pertur- 
bation from the hard-sphere model. A major in- 
convenience of using singular interactions such as the HS 
potential is that it is not possible to perform perturba- 
tion expansions in the small parameter pa, with p = N/L 
the particle density. Remarkably, the ground-state en- 
ergy of the strongly repulsive {~pa <C 1) LL gas coin- 
cides asymptotically with the ground state energy of the 
HS model, albeit with a < [3], that is, the eHS model 
with a < 0. With the novel momentum-dependent in- 
teractions developed in this work, we are able to perform 



perturbation theory for the strongly attractive or repul- 
sive LL gas, starting from the exact ground state for the 
eHS model, for both a > and a < 0. The success of 
such a perturbative treatment, which is seen a posteriori, 
is a priori expected, since the maximal difference of two- 
body T-matrices for the eHS and LL models at density p 
behaves as (Tohs — Ti^l)/ p = 0[(pa)'^], while the maximal 
difference between fermionized (Tonks-Girardeau) and LL 
T-matrices behaves as (Ttg — Ill)//0 = 0{pa), which is 
two orders worse in the small pa. These estimates show 
that, effectively, first-order perturbation theory from the 
eHS model corresponds to a third-order expansion from 
the fermionized Bose gas. 

The LL Hamiltonian ifLL can be written in terms of 
the eHS Hamiltonian i?cHS (Eq.- IU with momentum- 
dependent interactions ([5])) as 



LL 



oHS 



N 

E 

i<j=l 



S{xi - Xj)[gi^i, ~ g{kij)], (12) 



where gLL is the LL interaction strength and g(kij) is given 
by Eq. ([5]), with the limit - Xj — >■ 0+ as in Eq. ^ 
implicitly assumed. Making use of the bosonic symme- 
try of the particles, we can show that in any state of the 
eHS Hamiltonian, the expectation value E^^^^ of the per- 
turbation is given by -B^^^ = {S{xi - X2))J2i'<j=i9{kij), 
with {kij}ij the set of BA relative momenta of the par- 
ticular eHS (or HS) state, and g{kij) = gu^ — g(kij). 
In the thermodynamic limit (TL) i?^^) is well-defined 
for pa < 1/2, which sets an upper bound on the ra- 
dius of convergence of the perturbation expansion. We 
now apply the Hellmann-Feynman theorem, dE^is/da — 
J2i<j {^{xi~^j))dg{kij) / da, where ii'HS is the ground state 
energy of the HS Bose gas, given by [8] 



N' 



6m(l — pa)2 



(13) 



For the first-order energy correction we obtain i?'^^ 
C{p,a)dE}is/da, with 



C{p, a) 



(14) 



In the TL, the density of states (in the BA [3l]) for the HS 
Bose gas is a constant times a step function, and therefore 
we can write 



fldk^ p_dk2g{k,,) 



a{qa)^ 



" A *i A (fc-) " " 18 T(;a)2 ' (15) 



where q ~ 7rp/(l — ap), and where the approximate equal- 
ity is valid for q\a\ 1. 

The perturbative correction for the repulsive LL gas 
(a < 0) yields a minor improvement with respect to the 
eHS asymptotic value, so we concentrate on the attractive 
case, i.e. the sTG gas. In Fig. [l]we show the energy per 
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Fig. 1: Main figure: energy per particle as a function of pa 
for sTG gas in 1st order perturbation tlieory (black solid line), 
compared to HS result (dashed red line) and Monte Carlo re- 
sults of ref. [13] (blue dots). Inset: inverse compresibility for 
sTG to 1st order (black solid line) and for HS gas (red dashed 
line) . 



particle in the TL to first order in perturbation tlieory, 
E/N « (Shs + E^^^)/N , as a function of the gas parame- 
ter pa, and compare our results with existing Variational 
Monte Carlo (VMC) data from ref. [13] and with the HS 
result, Eq. (jl3p . Our results arc in excellent agreement 
with the VMC calculations up to pa ~ 0.25. Beyond that 
point and until pa w 0.3 our results deviate from VMC, 
but their tendency is still correct. For pa > 0.3 our cal- 
culation is not enough to describe the sTG gas, since in 
the present case the energy is overestimated in this re- 
gion. In Fig. [T] we also show the inverse compressibility 
mc^ = pdppL, with p the chemical potential, as a function 
of pa to first order in perturbation theory. Our results 
reproduce the overall features calculated in |13| and arc in 
good agreement until pa ~ 0.2, from where our calculation 
largely overestimates the fitted VMC results. 



Tan's contact for the sTG gas. — A quantity which 
has attracted much of recent theoretical and experimental 
interest is the so-called contact [H] , denoted by I. This is 
defined for bosonic and fermionic systems with zero-range 
interactions in any dimension pUll^l^BHSS] as the coeffi- 
cient of the asymptotic part of the momentum distribution 
no.(k), I = VL hmiki^oo fc^'T-cr(k), where 51 = is the vol- 
ume, D is the dimension and a is the spin component 
(omitted for spinlcss bosons). In the ID case, relevant 
here, it is related to short-distance correlations pUll32j . 
Relations between different properties of the system and 
the contact are generally known as Tan relations. Wc fo- 
cus here in Tan relations for the LL gas with or without a 
trap, which we then apply to the sTG gas. 

Tan relations can be proved in parallel to the higher- 
dimensional cases I1TJI35JI271I2H] by using the ID version 



of the so-called 77-selector pn[^[551 IM]. This reads 



ri[k) = 1 + 



-5{ll\k\). 



"i.9LL 

The energy of the system is given by 



(16) 



(17) 



with W = Y^^=i ^{xi) the total single-particle trapping 
potential. The adiabatic energy theorem reads 

dE _ h^I 
da m 

The two relations above were already known for the ho- 
mogeneous case [3JI20I. The generalized virial theorem, 
assuming W{x) 



(18) 



(X , is given by 



E 



13 + 2 



(W) 



— — c 

2m 



(19) 



Last, the pressure relation, which is only valid in the ho- 
mogeneous case, is given by 



PL = 2E 



h^Ia 
2m 



(20) 



In Fig. 121 we plot the contact per particle I/N for the 
homogeneous sTG gas obtained via perturbation theory 
from the eHS gas in the thermodynamic limit. As ex- 
pected, deviations of the perturbativc contact from the 
asymptotic HS result are more pronounced than for the 
energy as the gas parameter grows. The contact for a 
given quantum state of a system and, in particular, for 
a stationary state of the Schrodinger equation, is a pos- 
itive quantity. A violation of this property implies that 
the state under consideration is not a physical state of 
the system. As observed in Fig. |2l the contact exhibits 
a maximum at pa w 0.25, where perturbation theory is 
still approximately correct, showing that it will eventually 
become negative at a given critical value of the gas param- 
eter {pa)c where a super- Tonks-Girardeau gas cannot ex- 
ist. The VMC data of ref. [13] show that this is indeed the 
case, and from their fit one can estimate {pa)c ~ 0.4 — 0.45. 
For completeness, we note that, although our calculation 
is not quantitatively correct for pa w 0.25 — 0.3, the per- 
turbativc contact becomes negative at {pa)c ~ 0.38. 

In a realistic experiment, the sTG gas is created un- 
der harmonic confinement [14]. A qualitative picture of 
the trapped system can be inferred from the homogeneous 
Bose gas via the local density approximation (LDA) [35] . 
Within the LDA, we calculate the contact I for a harmoni- 
cally trapped sTG gas {W{x) = muj'^x'^/2), by making use 
of the virial theorem, Eq. ([TO)) . which is accessible with 
current experimental techniques |36j . We plot it in Fig. 
[2] as a function of iVa^/a^Q, where Cho = \/ h/moj. The 
appearance of a maximum for the contact and its subse- 
quent depletion - related to the gas-phase instability, as 
noted above for the homogeneous case - would be clear 
experimental signatures of the sTG gas. 
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Fig. 2: Main figure: contact per particle as a function of pa 
for sTG gas in 1st order perturbation theory (black solid line), 
compared to HS result (dashed red line). Inset: contact for 
harmonically trapped sTG to 1st order in the LDA from virial 
theorem, Eq. (|19|) . (black solid hne) compared to contact for 
trapped HS gas (red dashed line). 

Concluding remarks. — Wc have shown that the 
Bose gas with hard-sphere interactions is equivalent to 
a many-boson system with momentum-dependent con- 
tact intcrparticle potentials. The resulting model is 
Fourier-transformable, and constitutes, for the first time, 
a simple momentum representation for the singular one- 
dimensional hard-sphere potential. As an important ap- 
plication, we have used our equivalent, soft-core model 
as a starting point to obtain the properties of the attrac- 
tive and repulsive delta Bose gases in perturbation theory, 
and are in good agreement with the large-scale numeri- 
cal simulations of ref. |13| in the limit of applicability of 
our perturbation theory. Universal Tan relations for the 
Lieb-Linigcr gas are also derived and applied to estimate 
the contact for the super Tonks- Girardeau gas within our 
perturbation-theoretic approach. 

The momentum-dependent zero-range interactions in 
this work have been carefully defined and are general. 
Therefore, they are applicable to other one-dimensional 
systems: we can use them to construct integrable approx- 
imations to non-integrable systems from their two-body 
phase shifts, and treat the difference between the original 
and model interactions as a weak perturbation. 

* * * 
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